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We study quenched disorder in strongly correlated systems via holography, focus- 
ing on the thermodynamic effects of mild electric disorder. Disorder is introduced 
through a random potential which is assumed to self-average on macroscopic scales. 
Studying the flow of this distribution with energy scale leads us to develop a holo- 
graphic functional renormalization scheme. We test this scheme by computing ther- 
modynamic quantities and confirming that the Harris criterion for relevance, irrele- 
vance or marginality of quenched disorder holds. 
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I. IMPERFECTION 

There exist no perfectly ordered materials in nature. Every crystal's formation involves 
impurities sneaking in and getting stuck, lattice vacancies remaining unfilled: disorder is 
inevitable. These microscopic imperfections often leave conspicuous imprints on materials' 
macroscopic properties. A well-known example is Anderson localization in systems of non- 
interacting quasiparticles [1-3 J in which quenched disorder traps the quasiparticles, turning 
metal into insulating Fermi glass [I] . Similarly, quenched disorder in certain frustrated spin 
systems leads to glassy phases at low temperature [HI Ej. 

In strongly correlated systems, however, our theoretical understanding of quenched disor- 
der remains rather primitive. Is there a strongly correlated avatar of Anderson localization? 
What do Mott's law for direct current conductivity and its percolating picture in weakly 
correlated systems [7H9] morph into at strong coupling? Does many-body localization really 
happen [10J? What does quenched disorder do to systems governed by quantum critical 
points? Interesting theoretical questions abound. Meanwhile, from a pragmatic point of 
view, many technologically interesting systems, including the cuprate superconductors, are 
both strongly correlated and strongly disordered. It is clearly worthwhile to investigate the 
effects of quenched disorder in strongly interacting systems. 

To attack these challenging questions, we bring to bear holography, a powerful tool for 
studying the thermodynamic and transport properties of strongly correlated systems [TIV 
H3]. Our ultimate goal with this holographic approach is to find novel phases triggered 
by quenched disorder and to study transport properties within such phases. The goal of 
the present paper is more modest: we merely point toward promising places to look for 
interesting phenomena in the holographic context and begin developing some of the tools 
needed to explore them. Specifically, we trace both the flow of dilute disorder deformation 
and its effect on themodynamic quantities, working perturbatively around clean fixed points. 

As discussed below, self- averaging quenched disorder can be characterized by a distribu- 
tion Py [W(x)] over random functions W(x). Importantly, as we change the energy scale, the 
entire functional runs. In contrast to traditional setups where we have only a few relevant 
running parameters to keep track of, we must now deal with an uncountable infinity of run- 
ning couplings. To this end, we develop a holographic functional renormalization scheme, 
which enable us to compute disorder-averaged thermodynamic quantities in holographic 
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theories at finite temperature. 

We check the validity of our scheme by applying it to confirm that the Harris criterion 
around the clean fixed point holds for quenched disorder characterized by a Gaussian dis- 
tribution, at leading order in the strength of the disorder: the disorder is relevant when the 
scaling dimension of the strength of the disorder, u^is — d + 1 — 2Ae>, is positive, irrelevant 
when z/di S is negative, and marginal when this scaling dimension vanishes. Here, d is the 
spacetime dimension of the conformal field theory (CFT) and A© is the scaling dimension 
of the disordered operator. Whether marginal disorder is marginally relevant or irrelevant 
is a fascinating question we will revisit in a sequel [T5] by utilizing the technology developed 
herein. 

The organization of the paper is as follows. In SecjllJ we set up a holographic model with 
quenched disorder; our prototypical test case will be quenched electric disorder in CFTs 
dual to classical Einstein-Maxwell theory, but our results generalize straightforwardly. In 



Sec]III[ we study the bulk response to quenched-disordered boundary conditions, including 
backreaction on the metric, working perturbatively in the strength of the disorder. In 



Sec, IV, we then explore how the disorder distribution evolves as we change the energy 
scale. In SecjVj we propose a holographic functional renormalization scheme. Armed with 
the scheme stipulated in SecJVj we compute leading quenched-disorder correction to grand 



potential in Sec, VI In particular, we verify that the Harris criterion described above holds. 



We conclude in Sec, VII with a view towards the scenery beyond the perturbative regime. 



II. A MODEL HOLOGRAPHIC SYSTEM WITH QUENCHED DISORDER 

For ease of presentation, we will henceforth focus on a strongly correlated CFT which is 
holographic to classical Einstein-Maxwell theory with action 1 
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Here, the dimensionless constants 4? — 



^Gravity an d ^J— ' 



^Matter are determined by the 



parameters of the boundary CFT. We take a large ^Gravity limit to ensure classicality of the 



1 The special case with d = 1 + 1 should be treated with caution as the system would never enter hydro- 
dynamic regime. Our holographic calculations are performed for integers d > 2 + 1 and then results are 
analytically continued to all real numbers d > 1 + 1. 
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bulk theory, while keeping ^ Mattcr ~ 1 to bring the role of gravitational backreaction to the 

Gravity 

fore. The resulting classical equations of motion are 
1 



-d Q [V=gg QP g MN (d P A N - d N A P )] = (2.2) 



1 d(d-l) 8nG N 
ana jtCmn — mn ttto — 9mn 



FmpF]/ — -F P qF pq 'gMN 



(2.3) 



o"» M « or 2 yMJV n 2 

1 lL 9d+i 

where F MN = d M A N - d N A M . 

Let us now sprinkle impurities into a clean strongly correlated CFT defined holograph- 

ically as above. We will focus on the effects of a random "electric" potential V(x) = 

/ (2^-)d-i etk ' x ^(^) caused by quenched impurities in the system. In particular, the potential 

is time-independent: quenched impurities are, by definition, frozen on experimental time 

scales. The action of the clean boundary CFT, So, is thus deformed to 

S v = S + J dtd d - 1 xV(x)J'(f,x). (2.4) 

Here J^ 1 is a conserved U(l) current of the clean CFT which is dual to a bulk U(l) gauge 
field, A M . Via the holographic dictionary, this electric disorder translates into disordered 
boundary conditions on the bulk U(l) gauge field, 

lim A M (r,t,x) = 5 tM V{x). (2.5) 

Our choice of bulk coordinate system will be stipulated explicitly in the next section. 

Finally, the quenched random potential V(x) is assumed to self-average on macroscopic 

scales. 2 Given self- averaging disorder, we can legitimately estimate densities of extensive 

quantities, for example the grand potential Q, as 

" n ' 

in the thermodynamic limit where the volume of the sample Vd-\ approaches infinity. Here, 
Pv/[W / ( X )] is the functional associated with V(x) satisfying 

J VWP v [W]W(x 1 )...W(x n ) = ^-J ^- 1 x V(x 1 +x )...V(x n + x ) (2.7) 

and flw(x) is the grand potential of the system with an electric potential W(x). 



= fvWP v [W)(^) (2.6) 

d.a. J \ V d-1 J 



2 Physically, this means that homogeneity is approximately restored as we average measurements over 
regions much larger than typical disorder length scales. While the microscopic details of self- averaging 
disorder wash out, its effects persist via its statistical properties, for example in the rounding of the sharp 
Drude peak in real metals. To be sure, self-averaging is not universal, though it is very common. See, for 
example, Sec. III. A of [S] for a detailed discussion of when and why such disorder-averaged quantities give 
extremely accurate estimates of observable quantities for a macroscopic sample with quenched disorder. 
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A. Gaussian distribution and Harris criterion 

As a concrete example, let us consider the disorder dictated by a Gaussian distribution, 
randomly varying from site to site: 3 

JV[W(x)] = N,e-^tJ dd ^ W ^ = iV tt e-^k / (S^ W{k)W( " k) . (2.8) 

The normalization constant N$ ensures J VW Pv[W] = 1 and the constant /dis characterizes 
the strength of the quenched disorder. By dimensional analysis, /dis is seen to have the 
scaling dimension z/di S = d + 1 — 2A jt near the clean fixed point, where Ajt = d — 1 is the 
dimension of J 1 in the clean CFT. Thus, we would expect that the quenched electric disorder 
becomes irrelevant (relevant) at long distance if Ajt > (Ajt < %^), in other words, if 
d > 2 + 1 (d < 2 + 1). This Harris criterion arises, for example, in disorder- averaged vacuum 
two-point correlation functions. 4 We will use this criterion as a test on the machinery we 
develop in this paper to compute thermodynamic quantities. 



B. Power of holography 

It is in general not straightforward to compute disorder- averaged observables. For exam- 
ple, to obtain the disorder-averaged grand potential density, we must compute 

= J_ f VWP V [W} {In (Z w{x) ) } , (2.9) 

d.a. V d-l J 

which is not equal to y^Tn { / VW Py[W} [Z w ^) }. Computing the logarithm of the parti- 
tion function first and then disorder-averaging (not the other way around) is generally hard: 
dealing with this usually involves a handful of formal tricks, such as the replica trick and 
the cavity method. Looking through the holographic lens, our job is considerably simpli- 
fied by the fact that the logarithm of the partition function is automatically computed by 
the gravitational action of the disordered geometry. This, together with the holographic 
geometrization of the functional flow, makes holography a computationally tractable play- 
ground for studying certain aspects of quenched disorder in a strongly correlated CFT. 

3 This corresponds to uncorrelated impurities with J VW Py[W] W(x)W(y) — /di S £(x — y). 

4 An analysis for "classical" disorder has been carried out in [16] by using the replica trick. Interestingly, 
those results can be obtained by simply disorder-averaging without invoking replica at all. See also [T7j 
for holographic study of quenched disorder using the memory function method. 
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III. INTERMEZZO: BULK INFORMATION 



To systematically compute thermodynamic quantities and transport coefficients in this 
dirty holographic setup, we need to compute corrections to the geometry induced by the 
disordered boundary conditions. In this section we compute these corrections perturbatively 
in the strength of the disorder in two steps. First, we compute the leading bulk profile of the 
matter fields induced by their disordered boundary conditions. We then use this profile to 
compute a self-averaged matter energy-momentum tensor in the bulk, and use the resulting 
homogenous energy-momentum tensor to self-consistently compute the leading backreaction 
to the bulk metric at 0(fa s ). 

To set the stage, recall that the clean CFT at finite temperature and zero chemical poten- 
tial is described holographically by the following black brane geometry [18], here expressed 
in Euclidean time r = +it: 

g MN dx M dx N = +f(r)dr 2 + + ^ ^ dxlj (3.1) 

with f(r) = -[l--±y (3.2) 

The Euclidean time r has a periodicity so as to make the geometry regular. We 

must now compute the matter field profile subject to disordered boundary conditions in this 
undistorted geometry. 



A. Matter profile 



1. Maxwell dirt on pure anti-de Sitter (AdS) 



The quenched-disordered boundary conditions induce a nontrivial profile for the bulk 
£7(1) gauge field. Let us begin with the pure AdS solution where r + = 0. To first order in 



A M dx M 



(—idr). 



(3.3) 
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The bulk to boundary Green function Gq ( k; p ) is defined through 







, . (d-1) _k 



Go k;p 







v , (3-4) 
p-j \ / 

with G (k;p = ooj = 1 (3.5) 

and G (k;p = o) = 0. (3.6) 

The first equation is just the probe Maxwell equation of motion, the second is the asymptotic 



boundary condition (2.5), and the last is the requirement of regularity at the Poincare 



horizon. These can be exactly solved by 



Go k;p )=G 



r(¥) 



|k| 
P 



Ikl 

K d-2 

2 \ P 



where K d=2 (x) is the modified Bessel function of the second kind. 
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2. Maxwell dirt on hot black brane 



(3.7) 



For r + 7^ 0, we have 



A M dx M 



d 



■d-U 



ikx 



V(k) 



|g*bb ( 



L 2 k r 



(27i) d - 1 

with the bulk Green function Gbb ( k; p ) now satisfying 



r i r i 



(d 



G B b (k; p 

with Gbb (k; p = oo 
and Gbb ( k; p = 1 



-idr) 



1 

0. 



(3.8) 



(3.9) 

(3.10) 
(3.11) 



The last equality is the requirement of regularity at the black brane horizon. 5 For our 
purposes, the crucial property of Gbb ^k; pj , which will be key in showing the absence of 
inconsistent divergences in thermodynamic quantities below and which is derived in Ap- 
pendix |XJ is its high- momentum behavior near the boundary at p ^> 1: 

Gbb \ I „ / 1 



d p Un 



O 



for 1 < |k| < p. 



(3.12) 



Irregular solutions behave near the horizon as ~ I + K-(p — l)rn(p — 1) + 0(p — 1) 
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FIG. 1: An event horizon is distorted inhomogeneously by random electric field, giving rise to a 
rugged black brane. 

B. First-order backreaction 



The nontrivial probe profile of the gauge field introduces nontrivial energy-momentum 

tensor Tmn = ~4~ \FmpF n p — \FpnF PQ gMN\, seeding disorder nonlinearly into the right- 
s' l L J 

hand side of the Einstein equation. Its backreaction then reshapes the black brane geometry 
into the rugged one (see Fig.Q. Though it is a complicated task to obtain even the leading 
correction to the geometry for the random potential V(x), the algebra simplifies at long 
distance. Namely, as we zoom out to long distance in the x-direction, the quenched disorder 
self-averages and thus, in the bulk, an inhomogeneous energy- momentum tensor Tmn self- 
averages into the homogeneous one [Tj^jvld a ■ Similarly, an inhomogeneous rugged geometry 
9mn self-averages into the homogeneous geometry [<7Miv]d a a ^ ^ on § distance. 

To leading order, solving the Einstein equation with the homogeneous source [Tmn]^^ 
yields the resulting homogeneous [<?Miv]da • This is a straightforward, if stygian, exercise in 
gravitational perturbation theory, performed in Appendix [B] for the black brane geometry 



with the Gaussian disorder (2.8). In this paper, we need an explicit expression only for a 
rugged pure AdS solution with r + = 0, which we record here: 

f t 2 \ r_ „„„/rM , r, dr 2 



+ S){ i+ ^(i)}(8^ 



^r(i)} 



(3.13) 



,i=i 




FIG. 2: The flow of disorder can be represented by Feynman-Witten diagrams. Solid lines repre- 
sent probe propagations of U(l) gauge fields governed by Go in the vacuum whereas a wavy line 
represents a graviton propagation. 



with e = / dis x (|^) x (±)' 

xr\p) = 

and X2 Urc (p) = 



7i \ f (2d - 3)ai + (2d - 5)a 2 \ 

P ) 



d-1 

8tt 



d-1 J [(d-2 



(d-2)(2d-3) 
-2a 9 . I f p ,d-A 



dp p ' 



P-2 



(3.14) 
(3.15) 



Here, ai = / ^l^ y 2 {d y G^y)} 2 and a 2 = J ^ y 2 {Co(y)} 2 are constants of order 1, 
related to each other by (3d — 5)ai — (d + l)a 2 . 6 The integration constant p 2 will be 
appropriately chosen below. 



IV. FUNCTIONAL FLOWS 



The real utility of the holographic approach is that it geometrizes functional flows. Let us 
start with the pure AdS spacetime disordered by a random potential V(x) in the ultraviolet, 
characterized by a distribution Py[VT(x)]. To see how the distribution runs as we change 
the energy scale, we can evolve V(x) from infinity down to some hypersurface at r = r*. 
This gives us 14 (x) = zA,-(x, r*), from which we can read off the corresponding distribution 
Py^[W(x)] at energy scale ~ £|. 

Perturbatively, this process can be represented by Feynman-Witten diagrams (see Fig. [2]). 
In particular, holography provides an algorithmic way to keep track of functional flows. We 



6 To derive this relation, use Eq.(BlO) with r + = 0. 
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FIG. 3: Backreaction of disorder can also be represented by Feynman-Witten diagrams. 

also need to keep track of how the disorder distorts the pure AdS geometry at a given energy 
scale. This evolution, too, can be represented by Feynman-Witten diagrams (see Fig. [3]). 

A. Functional flow to first order in fa B 



Let us see how all this works at leading order in the strength of disorder. Formally flowing 
wn to r = r 
[see Fig. §»] 



down to r = r* convolves V(x) with the bulk-to-boundary Green function G {^-^ j , yielding 



= / (£^< k > G °(^) ' < 41 > 

In particular, for the case where the random potential is governed by the Gaussian distri- 



bution (2.8) at infinity, we obtain, at leading order in f^ s , 

PyJW(x)] = N^e' 1 ^ ^" (k,F( " k) with /*(k) = / dis |(? J . (4.2) 

We have already analyzed the first-order backreaction to the pure AdS geometry in the 
previous section, but there is one subtlety: choice of the integration constant p%. It specifies 
the scale at which we define "volume," which affects what we mean by "density." We will stick 
to the choice p 2 = oo for d < 2 + 1 and p 2 = for d > 2 + 1 so that e Xi (jj ^ /dis(yk) d_3 > 
choosing conventions around the ultraviolet and infrared stable fixed points, respectively. 
Similar remarks apply for "time" and "temperature," as usual. 

Note that, for d > 2 + 1, the disorder is irrelevant and thus starts to plague the ultraviolet 
geometry at large r. Thus we need to carefully renormalize disorder- averaged observables as 
we take the cutoff scale r c to infinity for this class of deformations. For d < 2+1, temperature 
provides an infrared cutoff scale, shielding us from infrared catastrophes. For d = 2 + 1, 
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while e Xi Ure (z) K /dfc, w e find a logarithmic running in [gu) d ^: toxT™ (z) 00 /disM^z;)- 
Whether this destroys the infrared geometry or the ultraviolet geometry depends on whether 
the deformation is marginally relevant or irrelevant. We will come back to this issue in the 
sequel [TJ] by going one order higher in f dis . 



V. HOLOGRAPHIC FUNCTIONAL RENORMALIZATION 

To compute disorder-averaged thermodynamic quantities at temperature T, we need to 
regulate them by cutting off the rugged black brane geometry at r = r c , then specify Dirichlet 
boundary conditions for the bulk fields - including the disorder functional - at the cutoff 
surface, evaluate the on-shell action, subtract temperature-independent divergences, and 
finally take the limit r c — > oo. Specifically, we propose a following recipe: 7 

1. At r = r c , fix the Dirichlet boundary conditions for the rugged black brane solution 
to be same as those of the rugged pure AdS solution with thermal-time periodicity 
Tperiod = z and volume V d „ x . 

2. Compute the regulated on-shell Euclidean action, which we identify as ^ip- via 
the holographic dictionary, for the rugged black brane solution. Renormalize away 
temperature-independent divergences by subtracting the r = r c surface contribution 
from the on-shell action for the rugged pure AdS solution. Take r c — > oo. 



A. Scheme in action to first order in /dis 



In the special case of quenched electric disorder governed by a Gaussian distribution (2.8) 
in the ultraviolet, this scheme works out as follows: 

1. We first find the rugged black brane solution whose 14(x) = (1 — -^H iv4 T (x, r c ) is 



7 For simplicity, we employ a background-subtraction scheme. A more systematic treatment as in [T5] would 
be interesting to pursue. 
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distributed according to [c.f. Eq.(4.2) 



PvJW(x)] = N lc e~ ! 2 /c« w W k ) W ith / c (k) = / dis G, 



L 2 k 



(5.1) 



and whose Euclidean time periodicity r per i d and volume Vd-i are set by 



''"period 



,,M 1 + ?^'°(r)U 



and / rf d x x 



d-1 



d-l 



1 + 



(d-l) 



-eoX 2 



pure / ' c 



L 



(5.2) 
^-(5.3) 



2. We then compute the on-shell action for the rugged black brane solution obtained 
above, subtract temperature-independent divergences, and take the limit r c — > oo. The 
Maxwell action contributes a temperature-independent divergence for the disorder- 
averaged grand potential density of the form 9 

Jd-ll, C / r2i \ ^ 2 



1 



9 2 /dis 



Gq 



while the contribution from Einstein action evaluates to 



d r Go 

Gn 



(5.4) 



(d- 1)L' 
SvrG 



r c ^d 



L- 



x 



^0 . pure / ' < 



i + fx? 



{d-l)e purc /rc 
2 X2 U 



(5.5) 



An easy way to get the latter is to note that, with the standard Gibbons- Hawking 
surface term, the only correction at 0(/dis) with respect to the undistorted pure AdS 
geometry comes from the change in the Dirichlet boundary condition at r = r c . 



On the field theory side, what the functional (5.1) at r = r c succinctly codifies is a com- 



plicated distribution entailing the cornucopia of multi-trace random deformations generated 
by integrating out the geometry from r = oo to r = r c [201 ETJ . Note that, due to the 



exponential decay of G^y) 



2 g 

r(^) 



y 2 Kd-2 (y) at large y, the momentum is cut off 



above Ikl ~ j% , conforming with the standard holographic intuition. 
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8 Note that we need to compare the random U(V) gauge field potential on the black brane geometry and 

the one on the pure AdS geometry in a properly rescaled time coordinate at r = r c , resulting in the extra 
i 

factor of ^1 — ^r^j above. 

9 Here, we are working with grandcanonical ensemble at zero average chemical potential. To work with 
canonical ensemble requires us to add an appropriate boundary term, changing Dirichlet boundary con- 
dition to Neumann boundary condition. 

10 Had there been no such effective cutoff, we would have suffered from infinite backreaction to the geometry. 
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VI. THERMODYNAMICS 



"CO 



d.a. 



We are now ready to calculate the disorder-averaged grand potential density 
to first order in f^, following the scheme developed above. The eventual satisfaction of the 
Harris criterion provides a nontrivial check on the machinery developed herein. 



Contribution from Maxwell 



Given the distribution (5.1), the 0(fd ls ) contribution from the Maxwell action to the 
disorder- averaged grand potential density becomes 

2 
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(lis 



d d ~ l k f fL 2 k 
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(6.1) 



Subtracting the temperature-independent divergence (5.4), we get 

2d- 3 
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d.a. 
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p d -'G % \ In ( ^) | - i-G 2 <9 P {In (G BB )} 



(6.2) 



p=Pc 



where p c = 

At low momentum, the integrand in the last parenthesis behaves well: for example, 



%) 



P=Pc 



approaches a finite number, (d — 2), since GbbI^- = 0; p 



1 — -3=2- However, for d > 2 + 1, contributions from high momentum modes with k ~ p c 
give rise to severe divergences: to see how the first term diverges, let us differentiate it with 
respect to p c : 



d" 



(2n) d - 1 
d d ~ l k 



p d - l G£d p { In 
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P Ctq 
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(6.3) 



where we used the property (3.12) in the last step. Similarly, the second term in the 
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parenthesis can be massaged into the form 
d^k 
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All in all, we end up with 11 
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,d-3 



(6.5) 



The divergent coefficient multiplying a temperature-dependent term looks horrifying at first. 
However, this divergence is precisely cancelled by a matching term in the Einstein action. 

B. Contribution from Einstein 



To evaluate the on-shell Einstein action for the rugged black brane geometry, we can 
employ the same trick we used for the rugged pure AdS geometry. The crucial step is to use 
the undistorted black brane geometry on-shell action 



(d-l)L d - 1 (T c \* r , 



8nG 



N 



(£) ^- x 1 



(d - 2) r d + \ 4tt L 2 



2(d- 1) r d j d r+ 



(6.6) 



with the modified volume V^'_ 1 = jl + ^^ eoXsT"" (^) } ^d-i an< ^ the horizon size r + set by 



4ttL 2 / H 



d r + 



(6.7) 



For d = 2 + 1, replace pzgyPc 3 by ln(p c ). 
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Subtracting the T- independent divergence (5.5) and taking the limit r c — > oo, we obtain 
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Upon using the identity (3d — 5)ai = (d+ l)a 2 , we see that the divergent term here precisely 



cancels the divergence we found in the Maxwell action, 



Vd-x 



d.a. 



C. Grand potential density: satisfaction of Harris 



With the naive divergences consistently cancelled out, we finally obtain 
-Q(TY 



d-l 



d.a. 



CO^Gravity^ + C l ^Matter /dis^ 3 + 0(/| is ) 



(6.9) 



where Co = — ^ (^) and Ci are constants of order 1. This formula happily accords with 
the Harris criterion. 



VII. IMPURE THOUGHTS 

In this paper, we have studied strongly correlated CFTs holographically dual to classical 
Einstein-Maxwell theory in the presence of the quenched electric disorder. In particular, 
we developed a holographic functional renormalization scheme and, for Gaussian disorder, 
computed the disorder-averaged grand potential density to first order in the strength of the 
disorder, /^is- The result accords with the Harris criterion, which adds confidence to the 
validity of our scheme. Namely, the quenched electric disorder dominates low-temperature 
thermodynamics for d < 2 + 1, whereas its effects essentially disappear at low temperature 
for d > 2 + 1. A particularly interesting case was d = 2 + 1 for which the quenched disorder 



was found to be marginal. In a forthcoming paper [15], we will investigate whether the 
quenched electric disorder is marginally relevant or irrelevant for d = 2 + 1. 
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These results indicate that holographic perturbation theory in /dis becomes unreliable at 
sufficiently low temperature for d < 2 + 1 (and possibly d — 2 + 1). Naturally, we would 
expect dramatic phenomena to emerge precisely when such perturbative analysis breaks 
down and the bulk geometry is significantly distorted. Thus we wish to embark on the 
journey beyond the perturbative regime, looking for transitions/crossovers lurking behind. 

Several powerful nonperturbative techniques have been developed in the study of spin 
glasses |5j |6]. For example, these techniques enable us to see the glass transition in the 
Sherrington-Kirkpatrick model for which classical mean field theory is valid. In particular, 
the replica method has been an extremely useful tool to analyze disordered systems in 
considerable generality. We wish to bring these techniques to bear on our particular problem 
at hand, which admits a dual representation in terms of classical gravitational theory. 12 

One natural scenario motivated by analogy to thermodynamic behavior of the 
Sherrington-Kirkpatrick model would be the following: for d < 2 + 1, as we decrease the 
temperature, the effective strength of the disorder fdi S T d ~ 3 grows, with perturbation theory 
breaking down when fdi S T d ~ 3 ~ 1. At this point, the system may enter into a glassy phase 
where we have not just one but many metastable solutions for a given asymptotic boundary 
condition V(x). It would also be exciting to see any connection to percolating picture of 
variable-range hopping for weakly correlated systems [9], for example a fragmentation of the 
black brane horizon. 
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Appendix A: Asymptotic expansion with Bessel envelope 



In this appendix, we analyze solutions to the probe Maxwell equation (3.9) in the black 
brane geometry in the high momentum limit. 13 To start with, we focus on the region very 
close to the boundary by going to the coordinate y 



M. 
p ' 



y y 



1 



1 _ vl_ 

|k| d 



(Al) 



If we neglect ^ for a moment, the equation is exactly solved by 



Giy) 



r(¥) 



y 2 Kd=2(y). 



(A2) 



For y >C 1, it approaches 1, whereas for y ^> 1, it exponentially decays as y 



2 e 



We will try to find a positive shrinking solution ^shrinking for Eq.( Al ) with s hrmking(k; y = 
0) = 1 which, for large |k|, rapidly decreases as we move away from the boundary at y — 0. 
There also exists a positive growing solution grO wing, say with the near boundary behavior 
</>growmg(k; ?/) = + 1 x y d ~ 2 + which rapidly increases away from the boundary for large 
|k|. The regular probe solution G B B w fth G BB ^k; p = ooj = 1 and Gbb(^', P = lj = 0isa 
linear combination of the two, but the coefficient in front of grO wing must be exponentially 
small for large |k| so as to satisfy the boundary condition at the horizon p = 1. Therefore, 
up to exponentially suppressed contribution, the high-momentum behavior of Gbb is entirely 
governed by (^shrinking- 

We now generate an asymptotic series for ^shrinking by using G as an envelope: 



-'shrinking (k, y) 

with ip (y) 



Giy) 
1. 



n=0 



(A3) 
(A4) 



13 For (and only for) d = 2 + 1, the result in this appendix can be reproduced by simpler WKB asymptotic 
expansion. 
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Plugging this into the probe equation yields a recursive relation 



/y yfd-3 
dy ' G${y') 

„.'d-3 ry' 



y 

- I dy' 

'0 



Go 2 (y') 



' n-l 



dy"GZ{y'y'nY,y" {n - l - l) W 



(A5) 



1=0 



Here the <i n 's are constants specifying the leading normalizable piece at each order in 
expansion. For a generic choice of d n 's, the corresponding solution grows exponentially for 
large y. Since we are looking for a shrinking solution, we will make a special choice of d n 's 
to tame such a rapid growth. Namely, we recursively choose 



^shrinking 



(d 



^ J™ dyG 2 (y)y 3 j^y^^ViCy) J • (A6) 



With this particular choice, we can inductively show that ip n (y) grows only as y( d+1 ) n 

for large y as opposed to generic exponential growth. In particular, the series provides a 

~ d 

nice asymptotic expansion as long as y <C |k| d +!. With the envelope, it follows that this 
special solution is in fact shrinking exponentially whereas generic solutions are exponentially 



growing. This shrinking solution (and hence Gbb) has the property advertised in Sec. Ill A 2 

I ^shrinking \ 



d p <!ln 



V Go J 



O 



p d-l k 2 



for 1 < |k| < p. 



(A7) 



Appendix B: First-order backreaction with the Gaussian distribution 

To first order in / dis , defining e + = / dis x \^-j2 ) x (l^Y 3 and 



we have 



**G N [T Tr ] d a . 
8nG N [T„] dA 
^G N [T^ 



in 



d- 1 



S7T 



d- 1 



d d ~ l k 



(2n) d 

d d - l \t 

(27i¥- 1 



ztKG 



BB 



|Cbb| 2 



(d- l)e. 

2L 2 
(d- l)e. 

2L 2 

(d-l)e. 
2L 2 



x f(r) x [-si(p) - s 2 (p)} , 
x y^y x [-si(p) + s 2 (p)] 



x — x 



(Bl) 
(B2) 



(B3) 
(B4) 

(B5) 
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with all the other components vanishing. As the self-averaged energy-momentum tensor is 
homogeneous, without loss of generality, we make the following ansatz for the self-averaged 
geometry: 



[g MN dx M dx N ] d . 



+f(r)U+e +X t 1 



BB 



dr 2 + 



dr 2 



+ 



L 2 



BB 



'd-1 



(B6) 



j=i 



We made a coordinate choice in r to set [g TT ] da = [5 ,rr ] da • Then, plugging it into Einstein 
equation and expanding to first order in e+, we get ordinary differential equations for x? B (p) 
and X2 B (p)- After lengthy manipulations, we arrive at following regular solutions: 



p 



dp' 



2(d- l)p 



/d-2 



d-2 



P 



12 



dp"s 3 (p" 



Pa 



and xTip) = ~ I dp' 



p> 



2 

7/2 



dp"s 3 (p") + 



{ Sl (p') - s,(p')} 



Pi) 



p'{(d-v-m^} 



(B7) 
(B8) 



with s 3 (p) 



p d {2(d - l)(2d - 3)p d -(d- 2)(d - 3)} si(p) 
{2{d-l)p*-{d-2)f 
p d {2(d - l)(2d - h)p d -{d- 2)(3d - 5)} s 2 (p) 
{2(d-l)p d -(d-2)} 2 



(B9) 



Here po is a free parameter related to a constant coordinate shift in r and p 2 is another 
integration constant. We can simplify the expressions further by using the identity 



1-4 



d 

dp 



pjMpy- s 2 (p)} 

2(d- I) 



(d-2) 



}J 



(BIO) 



which follows from Eq.(3.9) for Gbb(p)- Taking the limit r + — > of this solution yields the 



rugged pure AdS solution (3.13). 
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